Massive Dirac fermions and spin physics in an ultrathin film of topological insulator 
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We study transport and optical properties of the surface states which he in the bulk energy gap 
of a thin-film topological insulator. When the film thickness is comparable with the surface state 
decay length into the bulk, the tunneling between the top and bottom surfaces opens an energy 
gap and form two degenerate massive Dirac hyperbolas. Spin dependent physics emerges in the 
surface bands which are vastly different from the bulk behavior. These include the surface spin Hall 
effects, spin dependent orbital magnetic moment, and spin dependent optical transition selection 
rule which allows optical spin injection. We show a topological quantum phase transition where the 
Chern number of the surface bands changes when varying the thickness of the thin film. 
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I. INTRODUCTION 

A three-dimensional (3D) topological insulator is a 
novel quantum state of matter which possesses metal- 
lic surface states in the bulk energy gap.^"— The sur- 
face states consist of an odd number of massless Dirac 
cones which are protected by Z2 topological invari- 
ants. The first 3D topological insulator is Bi^^Sbi-^:, 
an alloy with complex structure of surface states which 
was confirmed by using angle-resolved photoemission 
spectroscopy (ARPES)i^"— Recently it was verified that 
Bi2Se3 and Bi2Te3 have only one Dirac cone near the 
r point by both experiments and the first principles 
calculation,—"— which attracts extensive attentions in 
condensed matter physics. The electrons or Dirac 
fermions in the surface states of topological insulator 
obey the 2-1-1 Dirac equations and reveal a lot of un- 
conventional properties such as the topological magne- 
toelectric effectJ^ It was also proposed that the surface 
states interfaced with a superconductor can form Ma- 
jorana fermions for performing fault-tolerant quantum 
computation^^"— Since the surface states surround the 
sample, it is still a great challenge for both experimen- 
talists and theorists to explore the transport properties 
for metallic surface states in the topological insulators. 

In this paper, we study an ultrathin film of topological 
insulator where tunneling between the surface states on 
the top and bottom surfaces opens a finite gap in the 
Dirac cone centered at the T point {k = 0). The low ly- 
ing physics of the ultrathin film can be described as two 
degenerate massive Dirac hyperbolas which form time 
reversal copy of each other. Each massive band has a k- 
dependent spin configuration: one near the T point deter- 
mined by the energy gap, and the other at k large enough 
determined by a spin-orbit coupling term quadratic in k. 
We show that the energy gap oscillates with the thin 
film thickness, and changes sign at critical thicknesses. 
Across the transition points, the k-dependent spin con- 
figuration near the F point is flipped while those at large 
k remains unchanged, leading to a topological quantum 



phase transition where the Chern numbers of the surface 
bands change [Eq. ((3T|) ] . In the two Dirac hyperbolas of 
opposite spin configurations, the k-dependent spin struc- 
ture results in a distribution of orbital magnetic moment 
and Berry curvature with opposite signs. In doped metal- 
lic regime, the Berry curvature drives the spin Hall effect 
of the extra carriers which leads to net spin accumula- 
tions on the thin film edges. We also discover a spin- 
dependent optical transition selection rule which allows 
optical injection of spins in the thin film. 

The paper is organized as follows. In Sec. |TT1 we 
present the derivation of the effective Hamiltonian for 
the thin film of topological insulator. In Sec. IIIIl the 
oscillation of the gap and the k-dependent spin config- 
uration are discussed. In Sec. IIVI the orbital magnetic 
moment and Berry curvature are addressed. In Sec. El 
the spin Hall conductance is derived in detail. In Sec. 
IVI[ the spin-dependent optical transition selection rule is 
shown. Finally, a conclusion is given in Sec. IVIII 



II. MODEL FOR AN ULTRATHIN FILM 

We start with the low-lying effective model for bulk 
Bi2Se3, in which surface states consist of a single Dirac 
cone at the F point<^ We take the periodic bound- 
ary conditions in the x-y plane such that k^ and ky 
are good quantum numbers, and denote the thickness 
of the thin film along z direction as L. In the basis 
of {|pl+,t),b2-,t),bl^,i),b2j,i)} which are the hy- 
bridized states of Se p orbital and Bi p orbital, with even 
and odd (— ) parities, the model Hamiltonian is given 

by 

where 
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and a a are the Pauli matrices, with k± = kx i iky and 
k^ = k^ + k^,. This model is invariant under time re- 
versal symmetry and inversion symmetry. In this pa- 
per, the model parameters are adopted from Ref. [H: 
M = 0.28cV, Ai = 2.2eVA, A2 = 4.1eVA, Bi = lOeVA^, 
B2 = 56.6eVA2, C = -0.0068eV, Di = l.SeVA^, 
£12 = 19.6eVA2. Though we adopt this concrete model 
to study the properties of an ultrathin film of topological 
insulator, the general conclusion in this paper should be 
applicable to other topological insulators. 

To establish an effective model for an ultrathin film, 
we first find the four solutions to the surface states of the 
model in Eq. ([1]) at the F point {k^ = ky = 0). 
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-A,) 



(2) 



where 
ho{A,) 



C ~ D^dl + (M + B^dl)a, ~ iAid.ax, 



The solution of the block-diagonal Hq can be found by 
putting a two-component trial solution into the eigen 
equation of the upper block 



hoiAi) 



^ E 



(3) 



with a, &, A the trial coefficients defining the behavior 
of the wavefunctions, and E the trial eigen energy. Note 
that the trial coefficients may have multiple solutions, the 
final solution should be a linear superposition of these 
solutions, with the superposition coefficients determined 
by boundary conditions. Then the problem becomes a 
straightforward calculation of the Schrodinger equation. 
Consider an open boundary condition that the wavefunc- 
tions vanish at the two surfaces located at z = zLL/2 of 
the film. Then we finally obtain two transcendental equa- 
tions. 



[C-M-E- {Di + Bi)Af]A2 _ tanh(AaL/2) 
[C - M - E -{Di+ Bi)\l]\i ~ tanh(Aai/2)' 



(4) 



note that where a =1 and 2,a=:2ifa = l, vice versa, 
so there are two transcendental equations. In Eq. (jlj, 
Ac define the behavior of the wavefunctions along z-axis, 
and are functions of the energy E 



K{E) 



(5) 



where for convenience we have defined F = Ai+2Di{E— 
C) - 2BiM and R = F^ - 4:{Df - Bf)[{E - Cf - M\ 
The self-consistent solution to the two equations in ^ 
can be found numerically, and give two energies at the F 
point, i.e., E^ and which define an energy gap 



A = 



E^. 



(6) 



Note that the bulk states with much higher or lower en- 
ergies due to the quantization of k^ in the finite quantum 



well in principle can also be obtained in this way, but are 
ignored here because we concentrate on the surface states 
near the gap. The eigen wavefunctions for i?+ and 
are, respectively. 



C4 



-(i5i + Bi>y+/+ 

-{D,+B,)ri-f_ 
lAJI 



(7) 



where C± are the normalization constants. The super- 
scripts of and 775*^2 stand for E± , and the subscripts of 
for parity, respectively. The expressions for and 
rif2 ^-^s given by 
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cosh(Aj''z) 



cosh(A^z) 



cosh(Afi/2) cosh(AfL/2)' 



sinh(Aj z) 



sinh(A2 z) 



sinh(A]^L/2) sinh(AjL/2)' 
(Af)2-(A±)2 



Af coth(Af i/2) - A^ coth(AjL/2) ' 



(8) 



where A^ = 
into Eq. © 



Ai tanh(Af L/2) - A^ tanh(AjL/2) 
\a{E±) can be found by putting back E± 

(b) 



T, = + 
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FIG. 1: (Color online) Schematic comparison between (a) 
the gapped K-K' valleys in the staggered graphene, and (b) 
the two-fold degenerate hyperbolas described by our effective 
Hamiltonian in Eq. (|10p . 



By replacing Ai by —Ai in the above solutions, the 
energies of the lower block h{){—Ai) of are found de- 
generate with those of /lo(Ai) and their wave functions 
are denoted as 'p^{—Ai) and x^{~^l)^ respectively. Now 
we have four states, namely, 0]"^, [x^(^i), 0]"^, 

[0, (/p^(-^i)]^, and [0,x^(-^i)]^, where f and i imply 
that the orbits with spin up and down are decoupled. By 
using these four solutions as basis states and rearranging 
their sequence following (note that each basis state is a 
four component vector) 











-A,) 











(9) 



we can map the original Hamiltonian to the Hilbert space 
spanned by these four states, and reach a new low-energy 
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effective Hamiltonian for the ultrathin film, 

Hcff — 



h+{k) 
/i_(fc) 



(10) 



in which 



^T^(fc) = Ef) — Dk^ — hvF{kx(Jy — kyffx) 
+T,{j ~ Bk^)a,. 



(11) 



Note that here the basis states of Pauh matrices stand for 
spin- up and spin-down states of real spin. In Eq. pil) . 
we have introduced a hyperbola index = ±1 (or ±). 
As shown in Fig. [TJ one can view the hyperbolas as the K 
and K' valleys in the staggered graphene [(a)], but being 
relocated to the T point [(b)]. Unlike the momentum cor- 
respondence in graphene, it is the CTz to — CTz correspon- 
dence in the present case. Therefore, the dispersions of 
h± are actually doubly degenerate, which is secured by 
time-reversal symmetry. Here, Tj = ± are used to distin- 
guish the two degenerate hyperbolas, h+{k) and h^{k) 
describe two sets of Dirac fermions, each show a pair of 
conduction and valence bands with the dispersions 



£,/„(k) = Eo~Dk^±y/{A/2 - Bk^)^ + {hvF)^k^, (12) 

where c and v correspond to the conduction and valence 
bands, respectively. The eigen states for Sc/v are 



"c/«(k) 



1 



-^c/v I 



(A/2 - i?fc2)r. 



(13) 



with \\u,,,\\ = ^[(A/2 - Bk-^)T, + £,/„J^ + {hvpy^k-^. 
Besides the gap A already defined in Eq. ([6]), the other 
parameters in Hamiltonian (|lip are given by 

D = (i32/2)((^t|a,|^^) + (x^|a,|x^))-i^2, 
B = {B,|2){{^^\a.^^^)-{x^\a^x^))^ 
Eo = {E++E^)/2, (14) 

and can be calculated numerically by using Eq. (O. 

The numerical results of A, vp, D, and B are pre- 
sented in Fig. [21 It is noted that \D\ must be less than 
otherwise the energy gap will disappear, and all dis- 
cussions in the following will not be valid. The A terms 
play a role of mass term in 2-t-l Dirac equations. 

In the large L limit. 



VF^{A2/h)^l-Dj/Bl 
The dispersion relation is given by 
£c/vik) = ±VFfik 



(15) 



(16) 



for small fc. As a result, the energy gap closes at fc = 0. 
The two massless Dirac cones are located near the top 



and bottom surfaces, respectively, as expected in a 3D 
topological insulator. 
In a small L limit. 



VF = A2/h, 



(17) 



and 



A = 2Bi7rVL2_ (|i8) 
The ratio of the velocity between the two limits is 



r, = l/Jl-Dl/Bl 



(19) 



It is noted that the velocity and energy gap for an ultra- 
thin film are enhanced for a thinner film. 
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FIG. 2; (Color online) [(a)-(c)]: Two-fold degenerate {tz = 
±1) energy spectra of surface states for thickness L — 20, 25, 
32 A(solid lines), and L ^ oo (dash lines). The grey area 
corresponds to the bulk states. The energy spectra are ob- 
tained by solving H{k, —idz)'i'{z) = E'i!{z) under the bound- 
ary conditions '^'{z = ±L/2) = 0. Please note that the scales 
of energy axis in (a)-(c) are different. The model parameters 
are adopted from Ref. i: M = 0.28eV, Ai = 2.2eVA, A2 = 
4.1eVA, Bi = 10eVA^ B2 = 56.6eVA^ C = -0.0068eV, 
Di = 1.3eVA^ D2 = 19.6eVA^ [(d)-(g)] The parameters for 



the new effective model Hcb: 
the Fermi velocity vf "vs L. 



D, B, the energy gap A, and 



III. ENERGY GAP AND K-DEPENDENT SPIN 
CONFIGURATION 

The opening of energy gap for the Dirac fermions is 
expected as a result of quantum tunneling between the 
surface states on the top and bottom surfaces. When 
the thickness of the ultrathin film is comparable with 
the decay length of the surface states into the bulk, the 
wavefunctions of the top and bottom surface states have 
a spatial overlap which leads to an energy gap at the F 
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point, analogous to the splitting of the bound and anti- 
bound orbitals in a double- well potential. The dispersion 
relation of the surface states are plotted in Fig. EJa), (b) 
and (c) for several thicknesses. A massless dispersion is 
obtained for a large L limit as expected. For a finite 
thickness, the energy gap at A; = is a function of the 
thickness L and decays quickly with L [see Fig. [IJf)]. It 
is noticed that the gap A even changes its sign at certain 
values of L. For instance, for the present case, at about 
integer times of 25 A. Correspondingly, the velocity of 
the Dirac fermions is also thickness dependent, which is 
enhanced for a small L. Strictly speaking, for an ultra- 
thin film, these so-called "surface states" emerge in the 
entire film. Nevertheless, they always lie in the bulk en- 
ergy gap and can thus be distinguished from those bulk 
originated quantum well states. 

From the solution to Eq. (|10l) . it is obvious that the 
spin vectors in each band take a k-dependent spin config- 
uration, in the neighborhood of the V point determined 
by the gap parameter A, and at large k determined by 
the term Bdzk^, as schematically illustrated in Fig. [31 
As the two Dirac hyperbolas are time reversal copy of 
each other, it is obvious that they have just the opposite 
k-dependent spin configurations, as can be seen from the 
Tz dependence of hr^{k) in Eq. pT|). 



(a) A / B < 



(b) A/S>0 



\ 



FIG. 3: (Color online) Schematic illustration of the k- 
dependent spin configurations in the conduction band (top) 
and valence band (bottom) in Dirac hyperbola = -f 1 when 
(a) A/B < and (b) A/B > 0. The center and corners of 
each panel correspond to that fc = and k is large enough, 
respectively. 



motion of the surface Bloch electron and is defined as 

e 



TO(k) 



E(k)\ |Vfeu(k)) -2,(20) 



where £(k) and u(k) are the dispersion and cigcnstates 
of hr,^>^ By putting Eqs. ([12]) and ^ into Eq. (pO]). 
we find the conduction and valence bands for each ft^^ 
have the same orbital angular momentum 



i(k) = 



h'^vj,{A/2 + Bk^) 



h 2[(A/2-Bfc2)2 



h^v^k^] 



(21) 



At the F point, the two degenerate Dirac hyperbolas ac- 
quire opposite orbital moments which add to the spin 
magnetic moment. 
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FIG. 4: (Color online) [(a)(b)]: Energy spectra of surface 
states (dash lines) and orbital angular moment m of the con- 
duction bands (solid lines) as functions of k for L = 20 and 
32 A. [(c) (d)] Berry curvature r2(fc) of the conduction band, 
(e) Orbital angular moment m of the conduction band as a 
function of thickness L for A;=0.01, 0.001, and AT^ . The 
singularities occur at where A changes sign. Only results for 
the hyperbola Tz = -|-1 are shown. The results for = — 1 
are right opposite to those for Tz = +1. 



IV. ORBITAL MAGNETIC MOMENT AND 
BERRY CURVATURES 

The opposite k-dependent spin configurations result 
in opposite physical properties of the surface Bloch elec- 
trons in the two Dirac hyperbolas, including the orbital 
magnetic moments and Berry curvatures, as shown in 
the context of graphenei^ These properties make pos- 
sible the manipulation of spin dynamics by electric and 
magnetic fields in the thin film topological insulator. 

Orbital magnetic moment arises from a self-rotating 



The k-dependent spin configuration also results in a 
gauge field in the crystal momentum space, known as 
the Berry curvature 



r!(k)-z-VkX (u(k)|iVfe|u(k)) 



(22) 



In an in-plane electric field, electrons acquire an anoma- 
lous transverse velocity proportional to the Berry cur- 
vature, giving rise to the Hall effect. In the conduction 
bands, the Berry curvature distribution near the F point 
is 



n{k) = -T, 



.(A/2 + Bfc2) 



2[(A/2-Bfc2)2+n2„|fc2]3/2 



(23) 
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The Berry curvature in the valence bands is right op- 
posite to that of in the conduction bands. There are 
opposite distributions of the Berry curvature in the two 
Dirac hyperbolas. The in-plane electric field can there- 
fore drive the spin up and spin down electrons towards 
the opposite transverse edges of the thin film, which is a 
surface spin Hall effect. 



V. SPIN HALL EFFECT AND TOPOLOGICAL 
QUANTUM PHASE TRANSITION 

The A term in our model plays a role as the magne- 
tization in the massive Dirac model exploited to study 
the anomalous Hall effecti^ In principle, we can find the 
Hall conductance for each h^.^ . Note that /i^-^ in Eq. pT|) 
can be explicitly written as 



Eq. (|27|) can be transformed into 



Eo - Dk^ 



i—x,y,z 



(24) 



where ai are the Pauli matrices, and the d(fc) vectors 
dx = hvpky, dy = -hvpkx, dz = '^zi^ ^ Bk^), 

(25) 

For the 2x2 Hamiltonian in terms of the d{k) vec- 
tors and Pauli matrices, the Kubo formula for the Hall 
conductance can be generally expressed as^i*^ 



e 

2h 



d'^'k {fk,c - fk,v) dda ddfj 



d3 



where d is the norm of {dx,dy,dz), fk,c/v — 
l/{exp[{ec/v{k) — fi)/kBT] + 1} the Fermi distribution 
function of the conduction (c) and valence (v) bands, 
with fj, the chemical potential, ks the Boltzmann con- 
stant, and T the temperature. 

At zero temperature and when the chemical potential 
lies between (— ■^), the Fermi functions reduce to 
fk,c = and fk,v = 1- By substituting Eq. (f25|) into ([26|) 
we arrive at 



g2 l-co 

'Jxy{o,Tz) = -T^^ d{k'^) 



{hvFYif + Bk') 
[{hvpYk^ + i^-Bk-^Y 



Note that by comparing above equation with Eq. ([23l) . 
we can write the Hall conductance in the form of the 
Berry curvature of the conduction band. 



0-0:2/(0, T^) : 
By defining 

cosO 



h J (27r)2 
(2 



i.fk,v - !k,c)^{k). 



^ BP) 



^{hv^Yk^ + ( 



(28) 



(29) 



'xy 



(0, 



Tz) = T, 

' 2h 



dik^) 



d{cos9) 
9(fc2) 



(30) 



The values of cos 6 at fc = and fc — > 00 only depend 
on the signs of B and A, respectively. As a result, in 
the insulating regime < M < we find that the 

anomalous Hall conductance for each hyperbola has the 
form 



CrxyiO,Tz) 



-^.^(sgn(A)+sgn(i3)), (31) 



where ^(sgn(A)-t-sgn(i3)) is the Chern number of the va- 
lence bands. It is a well known result in the field theory 
that the Hall conductance of the massive Dirac fermions 
is a half of quanta e^/h, (Jxy — — rz|^sgn(A) (if B = Q). 
(see Ref. [2J.) Our result demonstrates that a non-zero 
quadratic term in fc (i? 7^ 0) will give a reasonable result 
since a half quantized Hall conductance is not possible 
for a non-interacting system. For small L, the parameter 
B is always negative, and A changes its sign accompa- 
nying a gap close-and-reopen while the thickness of the 
thin film increases [Fig. ^f)]. The sign change of A 
fiips the k-dependent spin configuration near the F point 
and results in a jump Tz in the Chern number jS;26 j 
a topological quantum phase transition as discussed in 
quantum spin Hall effect in HgTe/CdTe quantum welk- 
in the ultrathin limit, A > and the Hall conductance 
is zero, while A < at a larger thickness and the Hall 
conductance ixyi^TTz) = Tz{e^ jh), as shown in Fig. [S] 
By the edge-bulk correspondence^^ gapless helical edge- 
states shall appear accompanying such a transition. This 
result is indeed supported by the solutions of the differ- 
ential equation of h±{k) in a geometry of a semi-infinite 
plane, in which there exists a gapless edge state only 
when A/B > 0, and the sign of the Hall conductance 
can be justified from the chirality of the edge states. 



L=20A (A/B<0) 



L=32A (A/B>0) 



(27) 
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FIG. 5: (Color online) The hyperbola-dependent Hall con- 
ductance vs Fermi level fj, for L — 20 A (a) and 32 A (b), 
respectively. 

When the Fermi level fi lies in the electron band, ac- 
cording to Eq. (US]) 
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where at zero temperature, 



Near the F point, the B term in the electron band dis- 
persion can be neglected, so that sa (^hvF)k^ + (^)^j 
and 



k p 



After a straightforward calculation and Taylor expansion 
of the result up to the linear term of /i, the Hall con- 
ductance is obtained when the Fermi surface is in the 
conduction bands, 



^.yip > y = -r, — [sgn(i?)+sgn(A) ] 



where the Fermi level fi is replaced by Sn, the carrier 
density near the band bottom, by using the relation 
|A|(/L(— ^)/{4:TTh'^Vp) Sn, which is found by using the 
dispersion Ec and the same Taylor expansion approach. 
One can check that when the Fermi level lies in the va- 
lence bands, axy is the same as the above result due to 
the particle-hole symmetry. 

As shown in Fig. Ele), B is negative for small L. For 
A > where the intrinsic film is a trivial insulator, dop- 
ing with electrons or holes can result in a metallic spin 
Hall effect, with the Hall conductance given by 



(32) 



where 5n is the carrier density. In the two Dirac hyper- 
bolas, the Hall conductance are opposite (see Fig. [5]) and 
spin vectors near the band bottom point in the -l-z and 
—z directions respectively. Thus net spin accumulations 
with out-of-plane polarization are expected on the two 
edges of the thin film. For A < 0, doping reduces the 
quantized Hall conductance and we find 



An 



(hvpfSn]. (33) 



The Hall conductance as a function of the Fermi level 
fi are plotted in Fig. [5] Note that because the Hall 
conductances for the two hyperbolas are always equal in 
magnitude and opposite in sign, and the two hyperbolas 
have right opposite spin orientations, here the Hall con- 
ductances are referred to as the ordinary and quantum 
spin Hall effects. 



the two gapped Dirac hyperbolas being time-reversal of 
each other, interband transition couples preferentially to 
right-handed (<7-|-) or left-handed (cr— ) circular polarized 
light, as was first discovered in the context of graphenej^ 
In the thin film topological insulator where the two Dirac 
hyperbolas are associated with opposite spin configura- 
tions (Fig. [3]), such an optical transition selection rule 
is of significance since it allows the spin dynamics to be 
addressed by optical means. The interband couplings in 
the two hyperbolas to normally incident circular polar- 
ized lights can be studied by calculating the degree of 
circular polarization, defined as 



k+(k)p-k-(k)p 
|^+(k)|2-f |^-(k)|2^ 



where ± corresponding to a± lights, respectively, and 
the interband matrix element of the velocity operator is 
defined by 

7r±(k)^(u,(k)|^±»|^K(k)), 

where Uc(k) and Mt,(k) are the eigen states for the con- 
duction and valence bands of h-r^ , and have already been 
given in Eq. (1131) . By ignoring Bk^ terms near the F 
point, one finds 

|7r±(fc)P~ (?ii;f)'(1±t,cos0)2, 

where cos6 ~ /S./{ec{k) — ev{k)), with Sc {£v) the dis- 
persion of the conduction (valence) band of hr^ given in 
Eq. (|12p . Then the degree of polarization is found out 
for each Tz 



V{k) 



2 cos 61 
1 + cos2 e ' 



Near the Dirac hyperbola center where £c{k) — e^ik) ~ 
|A|, a+ (a—) light couples only to the Dirac hyperbola 
Tz = sgn(A) {tz — — sgn(A)), causing transition between 
the spin down (up) valence state and the spin up (down) 
conduction state. This is in sharp contrast to the spin 
optical transition selection rule in conventional semicon- 
ductor which leaves the spin part of the wavefunction un- 
changed. This spin optical transition selection rule makes 
possible optical injection of spins in the thin film. For 
example, with band-edge optical excitation by a+ circu- 
larly polarized light, spin up electrons in the conduction 
band and spin up holes (by convention, a spin up hole 
refers to an empty spin down valence state) in the valence 
band are created in the Dirac hyperbola Tz = sgn(A), 
which can be separated by an in-plane electric field to 
prevent the radiative recombination. 



VI. SPIN OPTICAL SELECTION RULE 

Since an energy gap is opened in the Dirac hyperbo- 
las, interband transitions between the conduction and 
valence surface bands can be excited by optical field. In 



VII. CONCLUSIONS 

In this paper, we derive an effective model for an ul- 
trathin film of topological insulator, in which tunneling 
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between the surface states on the top and bottom sur- 
faces opens a finite gap in the Dirac cone centered at the 
r point (fc = 0). The low lying physics of the ultrathin 
film can be described as two degenerate massive Dirac 
hyperbolas which form time reversal copy of each other. 
Each massive band has a k-dependent spin configuration: 
one near the F point determined by the energy gap, and 
the other at k large enough determined by a spin-orbit 
coupling term quadratic in A:. It is found that the energy 
gap oscillates with the thin film thickness, and changes 
sign at critical thicknesses. Across the transition points, 
the k-dependent spin configuration near the T point is 
fiipped while those at large k remains unchanged, lead- 
ing to a topological quantum phase transition where the 
Chern numbers of the surface bands change [Eq. (|3T|)]. 
In the two Dirac hyperbolas of opposite spin configura- 



tions, the k-dependent spin structure results in a distri- 
bution of orbital magnetic moment and Berry curvature 
with opposite signs. In doped metallic regime, the Berry 
curvature drives the spin Hall effect of the extra carriers 
which leads to net spin accumulations on the thin film 
edges. We also discover a spin-dependent optical transi- 
tion selection rule which allows optical injection of spins 
by circular polarized lights into the thin film. 

This work was supported by the Research Grant Coun- 
cil of Hong Kong under Grant No.: HKU 7037/08P, and 
HKU lO/CRF/08. 

Note added: After posting this paper on arXiv, we 
learnt about the works by Liu et alM and Linder et alM, 
in which the similar finite size effect of the surface states 
was studied. 
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